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We investigate the possibility of preheating in hybrid inflation. This scenario involves at least two 
scalar fields, the inflaton field <j>, and the symmetry breaking field a. We found that the behavior 
of these fields after inflation, as well as the possibility of preheating (particle production due to 
parametric resonance), depends crucially on the ratio of the coupling constant A (self-interaction of 
the field a) to the coupling constant g 2 (interaction of 4> and a). For A>j 2 , the oscillations of the 
field a soon after inflation become very small, and all the energy is concentrated in the oscillating 
field 4>. For A ~ g 2 both fields a and <j> oscillate in a rather chaotic way, but eventually their motion 
stabilizes, and parametric resonance with production of x particles becomes possible. For A <C g 2 the 
oscillations of the field <j> soon after inflation become very small, and all the energy is concentrated 
in the oscillating field a. Preheating can be very efficient if the effective masses of the fields <j> and 
a are much greater than the Hubble constant at the end of inflation, since those fields can then 
oscillate many times per e-fold, with a large amplitude. Preheating can also be efficient if these 
fields are coupled to other light scalar (or vector) fields \- I n the recently proposed hybrid models 
with a second stage of inflation after the phase transition, both preheating and usual reheating are 
inefficient. Therefore for a very long time the universe remains in a state with vanishing pressure. 
As a result, density contrasts generated during the phase transition in these models can grow and 
collapse to form primordial black holes. Under certain conditions, most of the energy density after 
inflation will be stored in small black holes, which will later evaporate and reheat the universe. 
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I. INTRODUCTION 

In most inflationary models the energy density at the 
end of inflation is released when the scalar field oscillates 
near the minimum of its potential, creates elementary 
particles, which later collide, decay into other particles, 
and eventually thermalize via various interactions |Q. 
The original version of the theory of reheating of the 
universe was based on the idea that the oscillating scalar 
field can be considered as a collection of scalar particles, 
each of which decays independently ||. This simple 
picture, based on perturbation theory, is indeed valid 
in many inflationary models. However, recently it was 
shown that in many versions of chaotic inflation there is 
an additional stage when the oscillating scalar field pro- 
duces other fields in a regime of parametric resonance. 
This regime was called "preheating" ||. The process is 
nonperturbative and in certain cases it can be extremely 
efficient. It may change the final value of the tempera- 
ture of the universe after reheating [[M , it may lead to 
specific nonthermal phase transitions B, and it may pro- 
vide a new mechanism for the generation of the baryon 
asymmetry Q . For a detailed discussion of the theory of 
preheating see |t],|| and references therein. 

In this paper we will discuss preheating in the con- 
text of hybrid inflation models M . These models became 



rather popular lately because they can be relatively eas- 
ily formulated in the context of supersymmetric theories. 
They allow inflation to occur at very small values of scalar 
fields, and at a relatively small energy density, which re- 
moves some restrictions on inflationary models that pre- 
viously seemed almost unavoidable. This scenario is also 
interesting from the point of view of microwave back- 
ground anisotropies and large scale structure because it 
provides natural models for tilted primordial spectra of 
density perturbations |L0| . 

Hybrid inflation models describe at least two classical 
scalar fields, and a. After the end of inflation both of 
these fields oscillate, and the energy of their oscillations 
can be transferred to the production of (f> and a particles, 
or to the energy of other particles interacting with these 
fields. 

The investigation of reheating in these models should 
consist of several steps. First of all, one should study the 
oscillations of the fields (j> and a after inflation, neglecting 
particle production. Until now, this important question 
was only partially analysed in some special situations. 
As we will show, the behavior of these two fields after 
inflation crucially depends on the relation between cou- 
pling constants in the theory. In some cases soon after 
inflation one of these fields rapidly relaxes near the min- 
imum of the effective potential, whereas the other one 
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continues to oscillate. Then all the energy after inflation 
is concentrated in this one field, and reheating occurs due 
to its decay. We will show how to find which of the two 
field will play this role. In some other cases, the motion 
of these two fields becomes chaotic. Both of these fields 
oscillate simultaneously, and their energy transfers from 
one field to another many times during these oscillations. 
Without investigation of this issue one cannot achieve a 
proper understanding of reheating in hybrid inflation. 

The second issue to study is the possibility of a para- 
metric resonance which would transfer energy of the ho- 
mogeneous oscillating classical fields cj) and a to the en- 
ergy of elementary particles cf> and a. As we will see, this 
process in general is possible, but typically it is not very 
efficient. Then one should analyse a possibility of a res- 
onant decay of the fields (j> and a to some other scalar or 
vector particles. 

When this investigation is completed, one should study 
other channels of decay of the classical fields <j> and a, 
such as the usual (non-resonant) decay to their own 
quanta and into other particles in accordance to the ele- 
mentary theory of reheating Indeed, as it was shown 
in Ref. 0, this is a separate process, which differs from 
the resonant decay in the limit when parametric reso- 
nance becomes inefficient. 

Until now, investigation of reheating in hybrid inflation 
concentrated on the non-resonant decay of the fields cf> 
and a, see e.g. 11 . The results obtained by this method 



may be quite correct in some models where there is no 
parametric resonance, but if parametric resonance is pos- 
sible, then the theory of reheating becomes completely 
different. 

The purpose of this paper is to investigate this issue. 
Our paper is not going to be a complete study of pre- 
heating in hybrid inflation. We will restrict ourselves to 
the investigation of the first stages of preheating in hy- 
brid inflation, which is sufficient to identify the models 
where preheating may occur. A complete investigation 
of preheating in hybrid inflation can be achieved using 
the analytical methods developed in 0, together with 
numerical lattice simulations along the lines of Refs. |fl2| . 

The main part of our paper will be devoted to the 
standard version of hybrid inflation with one stage of in- 
flation. One may also consider inflationary models with 
two stages of inflation E3j. We will show that in such 
models both preheating and usual reheating are ineffi- 
cient. However, under some conditions reheating in such 
models may be rather efficient, but it occurs in a rather 
unusual way. In a recent paper |l4| ] we have shown that 
large density perturbations could be produced in such 
models at the time corresponding to the beginning of 
the second stage of inflation. In that case, large density 
contrasts re-enter the horizon during the radiation era 
and can collapse to form black holes. Radiation pressure 
may prevent the formation of primordial black holes, un- 
less their density contrast is of order one. However, if the 
process of reheating after inflation is very inefficient, then 
pressure after inflation will be negligible, and nothing will 



prevent the formation of primordial black holes. If these 
black holes have a mass of order 10 -5 — 10~ 9 g, they will 
evaporate very early and reheat the universe jljj ], pro- 
ducing all the entropy we observe today. This is a very 
interesting possibility which needs to be further explored. 



II. THE STAGE OF HYBRID INFLATION 

The simplest realization of chaotic hybrid inflation is 
provided by the potential B 
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The bare masses m and M of the scalar fields <f> and u are 
"dressed" by their mutual interaction. At large values of 
the fields, their effective masses squared are both positive 
and the potential has the symmetry a <-> — a. At small 
values of the field <p, the potential has a maximum at <f> = 
(7 = and a global minimum at <f> = 0, a = Co = M/y/X, 
where the above symmetry is broken. 

The complete equations of motion for the homogeneous 
fields are 



4> + 3H<j)= -(m 2 
a + 3H6 = (M 2 - g 2 



Xa 2 )a, 



subject to the Friedmann constraint 
H 2 



MI 2 



7 ~ l \ 2 + \a 2 + V{^a) 



(2) 
(3) 



(4) 



Motion starts at large <fi, where the effective mass of 
the a field is large and the field is sitting at the min- 
imum of the potential at a — 0. As the field cf> de- 
creases, its quantum fluctuations produce an almost scale 
invariant but slightly tilted spectrum of density pertur- 
bations 

During the slow-roll of the field <j>, the effective mass of 
the triggering field is m 2 = g 2 <f> 2 — M 2 . When the field 
4> acquires the critical value 4> c = M/g, fluctuations of 
the massless a field trigger the symmetry breaking phase 
transition that ends inflation. If the bare mass M of 
the a field is large compared with the rate of expansion 
H of the universe, the transition will be instantaneous 
and inflation will end abruptly, as in the original hybrid 
inflation model PJl7| . If on the contrary the bare mass M 
is of the order of H, then the transition will be very slow 
and there is a possibility of having yet a few more e-folds 
of inflation after the phase transition, see Refs. HI!. 

When cr = the inflaton potential becomes V(<fr) = 
M I AX + m 2 cf) 2 /2. Since the scalar field 4> is of the order 
of 4> c = M/g, for m 2 <C g 2 M 2 /X the energy density is 
dominated by the vacuum energy, 
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to very good accuracy Q. At that time a = 0, and one 
can neglect (f> in the equation of motion for the field 4>, 
so that 

3H <j> = -m 2 cj) . (6) 
It is then possible to integrate the evolution equation of 



(j)(N) = </> c exp(riV) 
,r 2 

r : 



m 



(7) 



where N = Ho(t c — t) is the number of e- folds to the 
phase transition. 

Quantum fluctuations of the inflaton field <j> produce 
long wavelength metric perturbations, 1Z = H6(p/<p, 
where 5<f> is the amplitude of the field fluctuation when 
it crosses outside the Hubble scale. These fluctuations 
give rise to a continuum spectrum of metric perturba- 
tions which can be computed exactly in the case of hybrid 
inflation lUl 



Vn = 



C(r) 2 g 2 M 2 _ 2rN 
r 2 6tt\M 2 



(8) 



where C(r) = T[3/2 - r]/2 r T[3/2], and the spectral tilt 
can be evaluated as 



n _ 1= d\nV n (k) = ^ 
din k 



(9) 



Note that the tilt is always greater than one in this model. 
Observations made on a wide range of scales, from COBE 
DMR to CAT experiments, impose strong constraints on 
the amplitude and tilt of the primordial spectrum J18| 



V]l 2 = 5x 1(T 5 (0.99 ±0.06), 
n = 0.91 ±0.10. 



(10) 
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In the limit m <C Hq, we have r ~ m 2 /3H 2 , see 
Eq. and C{r) ~ 1. This means that the amplitude of 
the curvature perturbation spectrum should satisfy fl6|| 



M 5 
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while the tilt of the spectrum is bounded by |lf 

7T M A 



(12) 



(13) 



III. THE END OF INFLATION AND THE ONSET 
OF THE STAGE OF OSCILLATIONS 

As we already mentioned, in the simplest version of the 
hybrid inflation scenario described above inflation ends 



as soon as the scalar field (j> decreases below <j> c — M/g. 
It will be important for us to investigate this process in 
a more detailed way, because it prepares the stage for 
the process of reheating which we are going to study. 
In particular, we should understand whether or not the 
"waterfall" process of symmetry breaking from a = to 

| = M/VX can be considered as a rolling down of a 
homogeneous field a. 

According to the classical equations of motion, the field 
(7 = cannot change its value because the first deriva- 
tive of the effective potential at a = vanishes. The 
process of spontaneous symmetry breaking in this case 
occurs due to the exponential growth of quantum fluc- 
tuations. Indeed, the field a has a (negative) effective 
mass squared — /i 2 (</>) = g 2 {4> 2 — </> 2 ), which vanishes at 
the critical point, but becomes large and grows up to 
/«(0) = M as the field <f> slides towards = 0. Quantum 
fluctuations of the scalar field a with momentum k grow 
as expwkt, where Wk = \f n 2 — k 2 and k — |k|. Sym- 
metry breaking occurs due to the growth of fluctuations 
with k < fi. This process produces an inhomogeneous 
distribution of the field a with (a) = 0. 

The resulting distribution of the scalar field a is rel- 
atively homogeneous on a scale / ~ /i -1 or even some- 
what greater |l9| because the rate of exponential growth 
is maximal at k = 0, and the fluctuations with k > fi 
do not grow at all. However, in our case we have 
an additional complication: the effective mass squared 



= g 2 



;) changes in time. In order to es- 



timate a typical scale on which the growth of the field 
a remains relatively homogeneous we will do the follow- 
ing. We will find /z 2 as a function of time At from the 
moment when <j) = <p c . The exponential growth becomes 
efficient when /x ~ At -1 . Later, fx(t) continues growing, 
and all modes o~k which began growing at the moment At, 
will continue growing exponentially with approximately 
equal speed, and eventually trigger a more rapid motion 
of the field (f>, see the next two sections. This motion 
will preserve initial homogeneity on a scale smaller than 
/i -1 ~ At, which we are now going to determine. 

The motion of the field <f> during this stage occurs 
due to the negative curvature of the effective potential, 
H 2 = g 2 {$ 2 c - <f) 2 ) » 2g 2 cj) c \<j>\At = 2M 2 m 2 At/3H = 
y / 2A/37rm 2 Afp At. Perturbations of the field begin to 
grow exponentially at At ~ /z , which gives 



y/Xm 2 M ¥ 



(14) 



Comparison of fx with H ~ M 2 /MpV~X shows that fj, >• 
H under the "waterfall condition" XmMp ^> Al 3 which 
should be satisfied in ordinary hybrid inflation ||. This 
is violated, however, in hybrid models with two stages of 
inflation 

This means that within the time [iT x <C H^ 1 the mass 
squared of the field a changes from being much greater 
than H 2 to being much smaller than — H 2 . Therefore 
there will be no specific inflationary fluctuations of the 
field a with the wavelength H^ 1 , but instead we will have 
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fluctuations of this field with the wavelength /i 1 , which 
will give us the typical scale of homogeneity I ~ yT 1 ~ 
(m 2 M P VX)- 1 / 3 . 

Now we should compare this length with other charac- 
teristic parameters which we will consider in our study of 
preheating. There are two orthogonal modes of oscilla- 
tions near the minimum of the effective potential, which 
we will encounter in our discussion. The first mode cor- 
responds to the oscillation of the field ct at cj) — 0; the 
corresponding frequency is given by m a = \/2M . One 
can easily verify that /x -C M for XMp ^> to 2 , which is 
satisfied in all hybrid inflation models. The only possible 
exception known to us is the model of Refs. HQ with 
two different stages of inflation; we will discuss reheating 
in this model separately. 

Another, orthogonal mode, corresponds to the fre- 
quency to.0 = gM/yX, see the next section. One can 
show that fj, < if to < lCT 3 A/pA 1/2 sr" 3/4 . This 
condition is also satisfied in most versions of the hybrid 
inflation model. 

The typical wavelength of perturbations produced by 
parametric resonance usually is not much different from 
the frequency of oscillations. Therefore our estimates 
indicate that in our investigation of preheating in the 
usual hybrid inflation one can treat the oscillating fields 
<f> and a as homogeneous. 



IV. PREHEATING 

Here we are interested in the behavior of the fields af- 
ter the end of inflation, as the fields oscillate around their 
minima, producing particles. The qualitative behavior is 
very different in two opposite limits, depending on the re- 
lation between the couplings A and g 2 , and independently 
of the relation among the bare masses, to <C H < M . 

We will show that explosive preheating of <p and ct par- 
ticles is not very efficient in hybrid inflation. However, 
we will consider an extra scalar field x, coupled to both 
d> and a, 



V(x) 



1 ,2 i2 2 i 1 ; 2 2 2 

2^1 ^ X + <* X 



(15) 



which will allow explosive production of x particles in 
hybrid inflation models for certain values of hi and h^- 

In order to study the production of particles of all fields 
we should analyze the equations of motion for the quan- 
tum fluctuations of the scalar fields. In linear perturba- 
tion theory we find 



8<t> k + 3H8'<l> k + (^ + «4) Hk = , 



a 2 
k 2 

2 

k 2 



5a k + 3HSa k + ^— + to 2 ) Sa k = , 



jik + 3H Xk + + m x) Xk = 0., 
where the effective masses are 



(16) 
(17) 
(18) 
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9 <* 

2i2 



1.2 2 



= 3 Act 2 + g z cp - M 2 + h\ X 
= h\i 



m 2 =h 2 4> 2 + hi ct 2 . 



(19) 
(20) 
(21) 



The rate of expansion of the universe is important only 
in the first few oscillations. However, it does play a role 
in ending the parametric resonance regime, since modes 
with physical momentum k/a in a resonance band will be 
redshifted by the expansion and fall out of the band, shut- 
ting off particle production in that mode. It is possible 
to take into account the rate of expansion by redefining 
the fields <p = a 3 ^ 2 S(j), tp = a 3 / 2 5a and X = a 3 / 2 x- Since 
the scale factor soon goes like a oc t 2 ^ 3 , the equations 
of motion for these fields are the same as Eqs. ( 16)-(p"8]) 
but without the friction term proportional to 3H. This 
will greatly simplify the analysis, while still taking into 
account the rate of expansion. 

As we mentioned in the previous section, there are 
two fundamental frequencies of oscillations near the min- 
imum of the effective potential: m a — \[2M and to^ = 
gM/y/X. Also, there are also two different scales of the 
fields 4> and a: 4> c — M/g, and ct = M/\f~X. 

If A ^> g 2 , one has cto <C 4>o and fh a ^> to^. In this 
case oscillations of the field ct tend to be insignificant. 
This field tend to follow adiabatically the position of the 
minimum of its instantaneous effective potential depend- 
ing on the value of the slowly changing field cj>. Most of 
the energy of these two fields will be concentrated in the 
oscillations of the field <f>. Thus the theory of reheating 
in this regime should describe decay of the field <j>. 

In the opposite limit, A -C g 2 , one has <jq ^> 4>q and 
to ct -c ?7i0. In this case the situation is reversed. Oscil- 
lations of the field 4> tend to be insignificant. This field 
tend to follow adiabatically the position of the minimum 
of its instantaneous effective potential depending on the 
value of the slowly changing field ct. Most of the energy 
of these two fields will be concentrated in the oscillations 
of the field ct, and the theory of reheating should describe 
decay of this field. 

The situation A ~ g 2 is more complicated. Both fields 
will oscillate with a comparable amplitude, transferring 
energy to each other in a rather chaotic way. 

The constants g and A could have in principle any 
value, as long as they satisfy the constraints ([l2]) and 
( |l3| ) put together, 
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M 



VA Mp 



(22) 



For M ~ 10~ 3 Aip, we have g 2 <C A, like in the usual hy- 
brid inflation ||, while for M <~ 10~ 16 Mp, we can have 
g 2 ^> A, like in some hybrid models with two stages of 
inflation |D|[l4|]. Finally there is a limiting case g 2 ~ A, 
which appears in the simplest models of hybrid inflation 
based on supergravity, see e.g. Refs. Jl7| , p0pl|Jll| ] . Re- 
heating in these cases is completely different from each 
other. We will analyze the three different cases in the 
following sections. 
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FIG. 1. The trajectory in field space (a, <f>) after the end of 
inflation, along the ellipse (4>/4>c) 2 + (°"/°"o) 2 = 1- Oscillations 
occur around <f> = and a = ao. This figure corresponds to 
A = l,g = 8 x I0~ 4 ,M = 10" 3 Mp,m = 1.5 x 10" 7 M P . 

A. Case g 2 < A 

There are two well differentiated regimes, depending 
on the ratio m^/H = gM/^/XH. For fh^/H 3> 1, the 
amplitude of oscillations of the inflaton field <f> after in- 
flation is large, of order one, while for fh^/H < 1, the 
amplitude could be very small. For both large and small 
ratios, the motion after the phase transition goes along 
the ellipse 

Act 2 + ,g 2 </> 2 = M 2 , (23) 

in the plane (a, <j>) , see Fig. [j]. As the field a grows to- 
wards its minimum at <ro, the field <p starts oscillating 
around zero, always along the ellipse. The main part of 
the oscillation comes from the (f> field (perpendicular to 
the a direction), at a ~ oq. 

4- = $(*) sin m,*t. (24) 

In the case of Fig. [I], equation (p3) is a good approx- 
imation for the oscillations of the inflaton field, where 
$0 — 1/15 and $(t) oc 1/t, see Fig. [| During the field's 
oscillations, the rate of expansion satisfies H(t) = 2/3t, 
which corresponds to a dust like era, ait) oc t 2 ^ 3 . 

One should note that the details of the behavior of the 
fields <p and <J may vary, depending on relations between 
different parameters of the theory. The results shown in 
Figs. |l| and || correspond to A = l,g = 8 X 1CT 4 ,M = 
10~ 3 M P - 10 16 GeV, to = 1.5 X lCT 7 Af P - 2 x 10 12 
GeV. In this case fh^ is of the same order as H at the 
beginning of oscillations. This is one of the main reasons 
why the amplitude of oscillations is so small. Note that 
the Hubble constant at the beginning of oscillations is 
approximately ten times smaller than Hq. It is a very 
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FIG. 2. The evolution after the end of inflation of H/Hq, 
4>/4>c and a /(To, as a function of the number of oscillations of 
the field <j>, N = m^t/Z-K. 

important circumstance, because it means that even be- 
fore reheating begins, the energy density drops down by 
about two orders of magnitude. There are essentially no 
oscillations of the field a, and the amplitude of oscilla- 
tions of the field (f> is relatively small. 

However, one may consider models with different pa- 
rameters, such that TO0 3> H . For example, one may 
consider a model with A = 1, g = 0.035, M - 10 12 GeV 
and to = 10 3 GeV. In this case ~ 10 5 H and the 
figures look somewhat different: The amplitude of oscil- 
lations of the field a becomes noticeable, but still remains 
very small, whereas the initial amplitude of oscillations 
of the field 4> becomes very large, comparable to <f> c . This 
allows for a long stage of production of <fi particles in the 
narrow resonance regime, before the backreaction and the 
expansion of the universe shifts the modes away from the 
resonance, see Ref. and below. 

It is possible to compute the production of particles of 
both fields during these oscillations. Let us start with 
the production of ^-particles. The effective mass ( |l9| ) of 
the field <j> along the trajectory (^) becomes 

2 2 , 9 2 M 2 ( 4> 2 \ , ,2 2 
m = m +h lX , 

~m 2 i (l-$ 2 sin 2 m^), (25) 

where we have neglected to 2 <C g 2 M 2 /X, see the discus- 
sion before Eq. (bh, and defined rh^ = gM/y/X. Initially 
(x 2 ) = 0, while the field (j> oscillates around zero, see 
Eq. (p^). Substituting into Eq. (|l6|), we arrive at the 
Mathicu equation for the production of <fi particles, 

y'l + [Af (k) - 2q cos 2z] <p k = , (26) 

where primes denote derivatives w.r.t. z = fh^t and 

k 2 

A^{k) = + l + 2q^ (27) 
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To investigate parametric resonance which leads to the 
production of particles one should find exponentially 
growing solutions of the Mathieu equation (|2^) in an 
expanding universe. A detailed explanation of this ap- 
proach to preheating can be found in Ref . [Q ; we will not 
repeat it here. We will only recall that preheating can be 
especially efficient in the regime when A^ik) < 2q§ and 
q<t> ^ 1/4, which corresponds to a broad parametric reso- 
nance. The resonance appears not for all k, but only for 
those k belonging to the instability bands of the Math- 
ieu equation. The concept of stability/instability bands 
is very useful if one can neglect the role of expansion 
of the universe; in a more general case one should solve 
Eq. ( p6| ) numerically and find out the modes which will 
grow specially fast. These modes will dominate particle 
production. 

For a complete investigation of preheating one also 
needs to investigate those effects related to backreac- 
tion of produced particles and their rescattering ]7],|],|l2| ■ 
We will make some comments on this issue later, but we 
will mainly concentrate on the first stages of preheating, 
where the corresponding effects can be neglected. This 
approach will allow us to identify those versions of hybrid 
inflation where preheating may be efficient. 

For the model described by Figs. ^ (with A = l,g = 
8 x fCT 4 ,M = 1(T 3 Af P ,m = 1.5 x 1CT 7 M P ) the initial 
value of q is extremely small, q^, = $q/4 ~ 10~ 3 <C 1. 
For k ~ 0, we are in the narrow region of the first res- 
onance band, A = 1. In this case there are solutions 
of Eq. (|2^) that correspond to exponential instabilities 
with occupation numbers of quantum fluctuations given 
by rik(t) oc exp(2/ifcz), which can be interpreted as par- 
ticle production ||. However, the effective parameter 
2fik ~ q will be extremely small and particle production 
will not be very efficient. 

It is possible to evaluate the time it takes for a par- 
ticular mode k to cross out of the resonance band, 
At ~ (ikH^ 1 , see Ref. ||. The occupation number then 
becomes nk ~ exp Explosive production 

shuts down when $ 4 < 32H/m$. In our model, this is 
true from the very beginning, and so we are confident 
that very little production of 0-particles will occur. 

In the models where the initial amplitude of the field 
4> is large (of the same order as <fr c ), production of <fr 
particles can be much more efficient. For example, in the 
model with A = 1, g = 0.035, M ~ 10 12 GeV and m = 
10 3 GeV one has case ~ 10 5 H, which means that the 
field can make about 10 4 oscillations before expansion of 
the universe becomes important. The initial value of the 
parameter q is about 1/4, see (psj). This corresponds 
to the narrow resonance regime in the second instability 
band of the Mathieu equation. 

The resonance in the second instability band at small 
q typically is very narrow and inefficient. Moreover, one 
should check carefully whether this band actually exists, 



because Eq. ( |26| ) is correct only when one neglects os- 
cillations of the field a. Meanwhile these oscillations are 
not entirely negligible, so they may affect the narrow res- 
onance in the second instability band. And indeed in our 
numerical investigation of the resonance in the coupled 
system of the fields <p an d a we did not find any evidence 
of 4> particle production. Even if this process happens, it 
is not really reheating; parametric resonance breaks the 
coherently oscillating field <f> into a collection of 4> parti- 
cles. A subsequent perturbative decay of <fi particles is 
necessary. 

Let us consider now the production of er-particles. The 
effective mass (pOj) along the trajectory (p3|) becomes 



ml = 2M 2 ( 1 - 



2 ,,2 



ml (1 - <f> 2 sin 2 m^t) , 



(29) 



where initially (y 2 ) = and m a = y/2M. Substituting 
back into Eq. (p7|), we arrive at the Mathieu equation for 
the production of a particles, 



< + [Afj (k) - 2q a cos 2z\ ip k = 0, 



where z = m^t and 



A a (k) 



(la 



9—2 1 
TO 2 $ 2 

ml 4 



-^f + 2q a 

mi. 



(30) 

(31) 
(32) 



The ratio m 2 /m 2 = 2A/# 2 , while q ~ (A/.g 2 ) $ 2 /2. Since 

<i> <C 2v / 2A/.9 in our case, we have that A — 2q ^> ^/q, i.e. 
far above the resonance band, and therefore we expect no 
explosive production of cr-particles from the oscillations 
of the field cf>. 

Moreover, a particles can neither be produced at the 
next (perturbative) stage of the decay of the field cf>. In- 
deed, the effective mass of ^-particles, ~ gM/*J\, in 
the case g 2 <C A, is much smaller than the mass of the a 
particles, \J~2M, so the decay of (f> particles is kinemati- 
cally impossible. Thus, in this regime the field cj> cannot 
transfer its energy to a particles. 

It is possible, however, to parametrically amplify fluc- 
tuations of another scalar (or vector) field x, which inter- 
acts with the fields <fi and a in accordance with Eq. . 
Since oscillations will occur around = 0, while a ~ <Jo, 
in order to have particle production in this field we re- 
quire that the induced mass from the symmetry breaking 
field er be much smaller than the corresponding oscilla- 
tions from the inflaton field </>, h\4> 2 ^S> h^cr 2 ,, which at 
the beginning of preheating corresponds to 



\h{ > g 2 h\ . 



(33) 



It is always possible to find couplings hi, hi that satisfy 
this constraint. Eventually, when the amplitude of oscil- 
lations of the <b field becomes 



G 



In n k 




FIG. 3. The exponential growth of the occupation number 
ilk of x-particles, for k — rh^, as a function of the number of 
oscillations of the field <j>, N = ifi^t/2-R. One can distinguish 
here the broad stochastic resonance regime. The typical value 
of the growth parameter is /i ~ 0.13 during the last stages of 
preheating. 



1*2 g 
hiVx' 



(34) 



the explosive production of x~P ar ticles will end. This 
may happen before or after backreaction sets in, depend- 
ing on the parameters of the model. 

The corresponding Mathieu equation for the produc- 
tion of x-particles is given by 



[A(k) - 2<jcos2z] X k = 0, 



where 



A(k) 



2q, 



h\M 2 $ 2 
g 2 frv^ 4 



(35) 

(36) 
(37) 



Initially, q = (Xhj/g 4 ) <5> 2 /A ^(Xhj/g 4 ) 10~ 3 , where we 
have assumed h 2 satisfies Eq. (|3^) and can be neglected. 
Since A 3> g 2 , it is possible to find natural values of 
the parameter hi that ensure go ^ 1- We are therefore 
in the broad parametric resonance region and explosive 
production of x-particles will occur. In this case there 
are solutions of Eq. (|3^) that correspond to exponential 
instabilities with occupation numbers of quantum fluctu- 
ations given by 



n k 



w-t( 



u] k (\X k \ 2 



\x k \ 2 



■ exp(2^ fc z) , 



(38) 



which can be interpreted as particle production ||] , where 
the frequency u k is given by 



2 j, 2 



(39) 



We have shown the exponential growth of the occupa- 
tion number n k of x-particles in Fig. [|a, together with 
the effective growth parameter fi k in the broad resonance 
region, see Eq. d38| ) and Fig. ||b. The model parameters 
used are A = = g = 8 x 10 _3 ,A/ = 10" 3 M P ,m = 
1.5 x 10 _7 Afp, for the typical momentum k — fh$. In 
this case, the typical growth parameter is relatively large, 
/j ~ 0.13, and particle production will be very efficient. 
Furthermore, for a large range of couplings we will enter 
the region of stochastic resonance, see Ref. [ItJ. 

The present case is analogous to the case found in 
for a massive inflaton field, where we substitute m — > 
fhfy = gM/y/X and g — * hi. All features found in that 
case will be present here, for a certain range of parame- 
ters. For instance, backreaction will set in at the time t\ 
given by M 



Ni 



2tt 



5 15 
8tt/i hi 



15, 



(40) 



when the expectation value of created particles satisfies 

,2 9 Ji /f 2 



g M 



(41) 



and the oscillating field 4> acquires a large effective mass. 
At that time the parameter q becomes 



1/2 

1/2 V 



hiV\ 
Ag 2 



10' 



(42) 



The efficiency of preheating at this stage is determined 
by the fraction of energy density in kinetic energy of x 
particles, ((Vx) 2 )i = fa^cj) 2 ^ g^ 2 , see Ref. Q, which 
could be very small compared with the energy density re- 
maining in the oscillations of the inflaton field, if q\ 3> 1. 

Even after backreaction there is still production of x~ 
particles in the broad resonance regime until the param- 
eter q falls below 1/4, see Ref. 0. At that time, 



( X 2 ) 2 ^ $2(t2) 



Ag 2 M 2 



■Qi 



1/2 



(43) 



where we have used the relation $2 ~ $1 q 1 between 
the amplitude of oscillations at both times, see 0. At 
this time the energy density is equally distributed be- 
tween the kinetic energy of the x particles and their inter- 
action energy with the field 0. In that case, ((Vx) 2 )2 = 



1/2 



see Ref. 



This result gets 



modified when including rescattering of x particles with 
particles of the field 4>. The end of preheating occurs 
slightly earlier, when $ r ~ 2.5<I>i q^ 1 ^ 4 , and the energy 
density is predominantly in the interaction energy be- 
tween x particles and the field (j), p x ~ \Q~ 2 g 2 cf) 4 ^ 4 . 

A relevant question in the presence of preheating is 
whether non-thermal effects due to large occupation 
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numbers of particles can restore the symmetry in the a 
field, see Ref. JtJ. For this to occur, we need that the 
effective mass in the false vacuum be positive at the end 
of preheating, 



V^'(O) = -M 2 



M 



hi (x 2 ) 



^g 2 h\ 1/2 
~Jh 2 ~ qi 



> o 



(44) 
(45) 



We see that the condition for symmetry restoration in 
this model requires 



( h 



1/2 



% /2 

— — > — ^s> 

4iVi Ag 2 h 2 



(46) 



see Eq. (|33|). Unless the initial value q is very large 
indeed, it is very difficult to satisfy this condition. There 
is however no fundamental reason why this should not 
happen, and, as a consequence, non-thermal symmetry 
restoration could in principle be possible in this model 
for certain relations between coupling constants. If this 
happens, the description of preheating given above will 
need to be modified. 

Let us now study the opposite limit, X g 2 , and see 
whether there is parametric resonance production of par- 
ticles in that case. 



B. Case A < g 2 

For A <C g 2 , the most important mode of oscillations 
is the oscillation of the field a near the minimum of the 
effective potential at er = M/\f\~. The initial stages 
of oscillations can be very complicated. Depending on 
the relations between the parameters of the theory, the 
first motion of the fields (<^>, a) may occur either in the 
(T-direction, as in the standard hybrid inflation scenario 
or in cj> direction, as in Ref. [jl3|,[l4|. However, in 
both cases the oscillations of the field cf> eventually be- 
come small, and all energy becomes concentrated in the 
oscillations of the field a, see Fig. 4. 

Thus, just as in the case A ^> g 2 , the motion of 
the fields eventually becomes essentially one-dimensional, 
but now it occurs in the a direction. The corresponding 
process is very similar to the one which occurs in new in- 
flation. The oscillations of the field a can be represented 
in the following way: 



1 + E(i) smm a t . 



(47) 



This is a good approximation for the oscillations of the 
a field with a small amplitude So <C 1, S(t) oc 1/t, see 
Fig. |. 

It is then possible to compute the production of parti- 
cles of both fields during the er-oscillations. Let us start 
with the production of cr-particles. Before the symmetry 
breaking field starts oscillating, the field has negligible 
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FIG. 4. The evolution after the end of inflation of H/Hq, 
and a /<jq, as a function of the number of oscillations 
of the field a, N — fh a t/2ir. This figure corresponds to 
A = 2 x 1Q- I6 ,g = 1.4 x 1(T 5 ,M = m= 1(T 8 A/ P . 



amp litude, so we will assume <fi ~ 0. The effective mass 
( pCf ) of the field a then becomes 



= 2M 2 + 3X(a 2 -a 2 ) + h 2 ( X 2 ) , 



(48) 



where (x 2 ) = initially, and rh 2 = 2M 2 . Substituting 
back into Eq. (|l7j), we arrive at the Mathieu equation for 
the production of a particles, 



^ + [A a (k) + 2q a sm2z] 
where z = rh a t/2 andQ 

2k 2 



0. 



(la 



2 M 2 



6E. 



(49) 

(50) 
(51) 



Initially in our model we have £ ~ 1/12, (x 2 ) = 0, so 
q a = 6 S ~ 0.5, see Fig. |. For k ~ 0, we are in the nar- 
row region of the second resonance band, A = 4. There 
will some production of cr-particles but very soon the 
modes will redshift away and production will shut down. 

The effective mass of (j> particles ( |l9| ) for small (f> and 
to <ti to^ becomes 



g 2 M 2 



hf (X 2 



771^(1 + 2E smm a t) , 



(52) 



where (x 2 ) — initially, and to^, = gM/^/X. This corre- 
sponds to the Mathieu parameters 



* We have neglected the S 2 and the $ 2 terms, which would 
have contributed as higher harmonics and soon be completely 
irrelevant. 
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2k 2 V 
a 2 M 2 + ~Y 



(53) 
(54) 



Therefore, for S < 1 one has -C A t f,(k). In this regime 
the resonance is very narrow, and one does not expect a 
large production of 4> particles. 

This conclusion remains valid for the perturbative de- 
cay of the oscillating scalar field a as well. Indeed, the 
effective mass of the field <f> particles ~ gM/\/~\ in 
the case g 2 ^> A is much greater than the mass of the 
a particles y/2M, so the decay of a particles to </> parti- 
cles is kincmatically impossible. Thus, in this regime the 
field <p cannot transfer its energy to a particles. This is 
very similar to what happens in the case g 2 <§; A, where 
the energy is concentrated in the oscillating field <fi, which 
(for g 2 <C A) cannot decay to a particles, see the previous 
section. 

Parametric resonance with production of x particles 
in principle could be possible. In the case where masses 
of particles near the minimum of the effective potential 
are not much greater than H, one can ignore the small </> 
oscillations. Then the \ particle production is described 
by the Mathieu equation 



Xk 



M 2 



hZcr, 



2h%a J E(t) sin fh a t] X k = , (55) 




A change of variables m„t = 2z 
I) to the Mathieu equation 



[A(k) - 2qcos2z]X k = 



7r/2 reduces 



(56) 



Here A(k) = 2k 2 /a 2 M 2 + 2h\/\, q = (2/if/A) E, z = 
fh a t/2, and primes denote differentiation with respect 
to z. 

Note that q < A(k)H. In this case the Mathieu plot of 
stability/instability bands shows that as soon as the am- 
plitude of oscillations of the field a becomes much smaller 
than <To, one has E « 1, which implies that q -C A(k). 
In this regime the resonance is narrow, and preheating is 
not very efficient ||. In our case E <C 1 from the very 
beginning of the period of oscillations, and x particle 
production is inefficient. 

It is possible in principle to conceive a scalar field x 
which couples not to a but to the difference a — erg, 



(57) 



In such models production of %-particle can be very effi- 
cient. Let us analyze this case is more detail here. 

Substituting @ into Eq. (||), we arrive at the Math- 
ieu equation for the production of x particles, (|35|), with 
z = That and 



A(k) 



h\M 2 
Xfh 2 



±2q, 

E^ 
4 ' 



(58) 
(59) 
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FIG. 5. The evolution after the end of inflation of 
H/Ho, (f)/4>c and a/ao, as a function of N = fha-t/2-K, in 
the case M ^> H. Note that the number of oscillations 
of the field a is about N/3, while the amplitude of oscil- 
lations of both fields is large. This figure corresponds to 
A = 1CT 2 , g = 1, m = 10 3 GeV, m = 1.3 X 10 11 GeV. 



Initially, ( X 2 ) = 0, so q = (h 2 /X)(E 2 /8) ~ (h 2 /X) lO" 3 , 
see Fig. ^. Since A <c g 2 in this case, we can easily get 
go > 1- This ensures explosive x _ particle production 
with a large effective growth parameter [ik ■ As an exam- 
ple, we considered a model with A = 2x 10" 16 , h — 10 g = 
1.4 x 10" 5 , M = m= 1CT 8 M P . We have found that the 
typical growth parameter is rather large, fi ~ 0.13, and 
particle production is very efficient. 

Note however that it is rather difficult (though not im- 
possible) to invent a realistic model where x particles 
have masses given by Eq. (p7|). We have found that it 
is possible in the context of supersymmetric models, but 
it requires explicit fine tuning of the superpotential of 
the theory. Thus, one may conclude that in a general 
class of hybrid inflation models where x particles acquire 
a mass /12O0 after spontaneous symmetry breaking, pre- 
heating for A <C g 2 is inefficient. This question should 
be reexamined in the theories with flat directions, and 
in the theories where x particles have small mass due 
to an accidental cancellation of two large contribution. 
Such a situation may not be entirely unrealistic; remem- 
ber that the Higgs doublet of the standard model has an 
extremely small mass as compared to the GUT scale. 

In the version of the hybrid inflation scenario which we 
discussed until now the amplitude of oscillations of the 
field a was relatively small. However, our conclusions 
concerning parametric resonance for A <C g 2 can change 
considerably if one finds the models with M 3> H , which 
induce a quick end of inflation and a large amplitude of 
oscillations of the field a after inflation. For example, one 
may consider a model with g = 1, A = 10 -2 , m = 10 3 
GeV, M = 1.3 x 10 11 GeV. This gives H - 2 x 10 4 GeV, 
?7i0 = 5 x 10 7 H and fh a = 8 x 10 6 H . In such a model the 
amplitude of the oscillations of the field a is very large, 
E ~ 1, see Fig. || Moreover, this amplitude will not 
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FIG. 6. The top panel shows the exponential growth of 
the occupation number riu of a particles, as a function of 
N — rhcrt/2-K, for k = 0.43M and the parameters of Fig. || 
It acquires a typical growth parameter /it, ~ 0.3 during the 
last stages of preheating. The lower panel shows the occu- 
pation number nk of <j> particles, for k = 0.2M. The growth 
parameter is an order of magnitude smaller, Hk — 0.023, in 
this case. The dashed line shows the oscillations of the field 
a. The upper panel shows lOcr/cro, while the lower panel 
shows a /ao. The number of the a field oscillations differs 
from TV = fh a t/2ir because the oscillations are not harmonic. 

be damped by the expansion of the universe during 10 6 
oscillations of this field. As a result, in this model there 
will be a very efficient production of a, 6 and \ particles 
in a broad resonance regime, see Figs. ^ and 0. 

These figures illustrate the growth of fluctuation ne- 
glecting backreaction. In fact, preheating in this model 
is so efficient that it completes within few oscillations. 
Indeed, let us estimate how strong should be the growth 
of rik in order to transfer all energy of the oscillating 
scalar field a to a particles (this is the leading process). 
The initial energy density of the field a is M 4 /4A. Each 
produced particle carries the energy ~ fh a ~ V2M. This 
implies that one should have n k ~ A -1 in order to trans- 
fer all energy of oscillations into a particles. The total 
duration of the process in terms of the number of oscilla- 
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FIG. 7. The top panel shows the occupation num- 
ber nk of x particles as a function of N = m a t/2n, for 
k — 0.2M, hi — 0.3, hi = 1 and the parameters of Fig. ^. 
It acquires a typical growth parameter fx ~ 0.03 during the 
last stages of preheating. The lower panel shows the same as 
above but for hi = hi = 1. The growth parameter is some- 
what smaller, fj, ~ 0.013, in this case. The dashed line shows 

0-/(70. 

tion N follows from the relation exp(47r/iiV) ~ A -1 , which 
gives N ~ (47r/i) _1 In A -1 . For (j, ~ 0.3 and A — 10~ 2 one 
has N = O(l), i.e. the whole process finishes within one 
or two oscillations; the part of Figs. || and [7] at TV ^> 1 is 
redundant. 

The unusual efficiency of the process of a particle pro- 
duction can be explained by the relatively large value of 
the coupling constant A and by the tachyonic instabil- 
ity of the field at the first stages of its rolling down. A 
detailed theory of this process will be considered else- 
where |jg|| . A simple interpretation of our results is as 
follows. 

The process of falling down of the field a in the regime 
A <C g 2 occurs almost independently of the field 4>. The 
field (j> triggers the process, but then the field a moves 
essentially as a field falling down from the point a = in 
the theory with the effective potential ^(u 2 — ctq) 2 - As we 
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discussed in Section III , neglecting expansion of the uni- 
verse, fluctuations of the field a in the beginning of the 
process grow as exp \Jk 2 — M 2 t. Initially, there was no 
homogeneous classical field a, the whole process of spon- 
taneous symmetry breaking occurs due to the growth of 
quantum fluctuations. The value of the homogeneous 
mode (a) vanishes all the time; spontaneous symmetry 
breaking occurs when yj {a 2 ) created by the instability 
approaches cr . As a re sult, already after the quarter 
of the first oscillation, y/ (a 2 ) approaches <Jo, all energy 
will be concentrated not in a homogeneous field a but in 
its fluctuations with a wavelength I ~ 2nM~ 1 In 1 / 2 
which is somewhat greater than M _1 113]. In this sense 
one may even say that preheating ends within one quar- 
ter of an oscillation. This, however, would be not quite 
correct because after rolling there will be a stage of res- 
onant production of particles with momenta comparable 
to M. 

We should emphasize, however, that here again we deal 
not with the decay of the field a to other particles, but 
with a rapid process of transformation of the energy of 
the classical field a to the energy of a particles. The 
subsequent decay of these particles occurs in accordance 
with the elementary theory of reheating ||. 

Let us now briefly discuss the case when the two cou- 
plings g 2 and A are of similar magnitude. 



C. Case g 2 ~ A 

This case corresponds to the simplest supergravity hy- 
brid inflation model Jl7],^o],^l] , based on a superpotential 
of the type 



W = S(k< 



(60) 



which induces a scalar potential with self-coupling A = k 2 
for the triggering field 0, and couples the two fields with 
strength g 2 — k 2 , see Eq. (0). Thus, we have the same 
coupling, A = g 2 , as one should expect from supersym- 
metry universality.]]] 

In this case, the motion is always two-dimensional in 
the (a, 4>) field space: the field <f> oscillates around zero 
while the field a oscillates around <7q , with similar ampli- 
tudes but varying frequencies, see Fig. |^. Since the fields 
are coupled, their effective masses also vary. This in- 
duces the chaotic behavior of the fields shown in Fig. 
We have taken as parameters, A = g 2 — 10~ 6 ,Af = 
1(T 6 M P , m = 2 x 1CT 10 M P . The amplitude of the fields 



^This is not a general rule. For example, it is not necessarily 
true in D-term hybrid inflation S|, where the couplings g 2 
and A come from the F and the D terms respectively. 

'For a general discussion of chaos in hybrid inflation and in 
preheating, respectively, see Refs. [^3) and jjj,^4|. 




H 




1.5 



0.5 



-0.5 



N 

FIG. 8. The top panel shows the motion of the two scalar 
fields in the plane (a/ao, <j>/<f> ). The lower panel shows the 
evolution after the end of inflation of H/Hq, <f>/<f> c and a/ao, 
as a function of the number of oscillations of the field a, 
N — fhcrt/2-K. It is clear that for the first dozens of oscil- 
lations, the behavior of both fields is extremely chaotic, and 
the amplitude of the fields decreases rather slowly compared 
to previous cases. These figures correspond to the parameters 



A 



10" 6 ,M = 10"° M P ,m = 2x 10" 



1 M P 



takes many oscillations to decrease significantly, e.g. af- 
ter N = 50, the amplitude of the <j> field is still of or- 
der 1/2, although after one e-fold (corresponding here 
to N ~ 340 oscillations of the a field) it decays like 1/t 
with amplitude <f>o ~ 1/5- Meanwhile its frequency is 
rather chaotic for the first e-fold, as one can appreciate 
in Fig. |^, but eventually regularizes and becomes oscil- 
latory, while the rate of expansion approaches the usual 
dust dominated behavior, e = —H/H 2 = 3/2. This can 
be appreciated in Fig. |^. 

Even though initially the chaotic behavior of the fields' 
oscillations prevents any particle production, as soon as 
their behavior becomes regular (after about 1 e-fold) the 
amplitude of the oscillations is large enough (thanks to 
the fact that it has not decreased significantly during 
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120 



FIG. 9. The evolution of the parameter e = —H/H 2 after 
the phase transition in the case g 2 — A. The dust dominated 
era, e = 3/2, is approached rather chaotically after many field 
oscillations, although in less than one e-fold. 



the chaotic motion) that one may have some parametric 
resonance. 

Like in the previous case (g 2 S> A), a new field x could 
couple to a and <fi with potential (|l5|). Since a acquires 
the vacuum expectation value uq , we have to suppress the 
coupling to a in order to induce parametric resonance of 
X-particles. In this case, since g 2 — A, the condition 
( p3| ) becomes hi ^> hi. Eventually the amplitude of the 
field 4> oscillations will decrease to $(i) < h-xjh\ and the 
explosive production will end. The Mathieu eq uation 
is the same as in the previous case, see Eq. (p5[), with 
occupation number ( |38| ) and frequency (|39|). 

We have taken in this case hi ~ and hi = 4 x 10~ 3 , 
and found explosive production of x-particles, see Fig. |lC| . 
Note the initial plateau, as the chaotic behavior of <j> pre- 
vents a significant production of x-particles. Eventually, 
the stochastic resonance regime begins, but the growth 
index jik — 0.03 associated with the mode k = M at the 
last stages of reheating is not as large as in the previous 
cases. 

Backreaction of the produced x-particles on the oscil- 
lations of the </>-field will become important at time t\, 
given by Eq. (Eoj), when the expectation value of the x~ 
field becomes 



U h ~ h\ h 2 ' 



(61) 



and the field cf> acquires a large mass ~ h\ (x 2 )i- Particle 
production still continues until the effective q parameter 
enters the narrow resonance regime, q <■ 1/4 Q|. 

Symmetry restoration, if it exists, depends only on the 
couplings of x to and a, see Eq. (pi|), 



(62) 
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FIG. 10. The exponential growth of the occupation num- 
ber rik of ^-particles as a function of the number of oscillations 
of the field a, N — fh a t/2n, for k — M and the parameters 
of Fig. [| One can distinguish here the chaotic plateau fol- 
lowed by the broad stochastic resonance regime. It acquires 
a typical growth parameter \x ~ 0.03 during the last stages of 
preheating. 




N 

FIG. 11. The exponential growth of the occupation num- 
ber rifc of a-particles as a function of the number of oscillations 
of the field a, N = m CT t/27r, for k = M and the parameters of 
Fig. |^. One can only distinguish here the narrow resonance 
regime. It acquires a typical growth parameter jj, ~ 0.01 dur- 
ing the last stages of preheating. 



In this special case, there is also a small production of 
cr-particles, due to the oscillations of the homogeneous a 
and (/(-fields. The effective mass ( pfj|) of the cr-fluctuations 
in this case is 



2 ^ 2 + 3\(<j 2 -a 2 ), 



m^+g 
m 2 + M 2 ® 2 sin 2 fh^t 

+ 3M 2 ( 2S sin m a t + S 2 sin 2 fh a t 



(63) 



(64) 



Although difficult, this condition could still be satisfied. 



where we have used Eqs. ( |24| ) and ( ^7| ) 

Since m| = g 2 M 2 /\ ~ M 2 and m 2 = 2M 2 , the ef- 
fective mass (p3[) oscillates with various frequencies. In 
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this case, we cannot neglect the terms proportional to 
<fr and £ 2 , since they are not small and represent the 
oscillations of the two background fields. Initially, the 
motion is chaotic and no parametric resonance is possi- 
ble. Eventually, oscillations become more regular, when 
the amplitude of oscillations is still large, £ ~ 1/2 and 
$ ~ 1/3. Then the Mathieu equation for the production 
of a particles becomes 



ip'l + K(fc) + 2q a sin 2z\ ip k = 0, 



where z = rria-t/2 and 
A a {k) -- 
q* -- 



4fc 2 



12M J 



+ 4, 



(65) 

(66) 
(67) 



Initially we can neglect backreaction, and m 2 = 2M 2 . In 
this case q a = 6Eo ~ 3, see Fig. ||d. There will be a slight 
production of c-particles, as can be seen in Fig. Oa. 
However, the effective growth parameter /i^ ~ 0.01 is 
never large enough to dominate the decay of the er-field. 
Backreaction will be dominated by the x _ P ar ticle produc- 
tion, certainly not by cr-particle production. 

The existence of the regime of chaotic oscillations may 
be important in application to realistic models of hy- 
brid inflation. For example, recently it was argued that 
in hybrid inflation scenario based on the SU (5) model, 
spontaneous symmetry breaking has a tendency to oc- 
cur in the direction with a "wrong" type of symmetry 
breaking SU (4) x U (1), which was considered as a prob- 
lem of this scenario |25j| . However, analogous result is 
true even for the standard high temperature phase tran- 
sitions in the SU(5) [^6), where the phase transition 
is typically first order, and the bubbles of the phase 
SU(4) x U(l) form much more often than the bubbles 
of the phase SU(3) x SU{2) x U(l) g7|- Still, the uni- 
verse may eventually end up in the proper vacuum state 
SU{3) x St/ (2) x U(l), after the bubbles of the phase 
SU(3) x SU{2) x U(l) "eat" the bubbles of all other 
phases J2?],[l| . Similarly, the phase transition in the SU (5) 
version of the hybrid inflation model may be first order, 
with several different phases being produced. Then one 
should study the equilibrium between these phases. But 
even if the phase transition is second order, as in the sim- 
plest model considered above, one may wonder whether 
the chaotic behavior of the scalar fields <f> and a may split 
the universe into many domains with different types of 
symmetry breaking. Then one should study which of 
these domains will eventually survive. This process be- 
comes even more complicated if one studies backreaction 
of x particles produced during preheating on the shape 
of the effective potential. The best way to investigate 
this issue is to use lattice simulations Jl2],|l!| . 



V. REHEATING IN HYBRID MODELS WITH 
TWO STAGES OF INFLATION 

So far we have studied the generic behavior of pre- 
heating in various hybrid inflation models. We see that 
in some of them it is indeed possible to produce large 
amounts of particles in an explosive way and thus re- 
heat the universe rather efficiently. In all models which 
we studied so far inflation ends at the moment when the 
field 4> reaches its critical value 4> c and a phase transition 
with the generation of the classical field a occurs. 

However, this is not the most general situation. There 
are some hybrid inflation models where inflation ends be- 
fore the phase transition, see e.g. Ref. (2^^l[]. There 
is also an interesting class of models where inflation 
continues for a while after the phase transition [ |l3| . |l4[ . 
This happens if the symmetry breaking field has a vac- 
uum expectation value Co ~ Mj>, and a mass of order 
M ~ 1 TeV. This requires an extremely small parameter 
A <C g . As a consequence, after the phase transition at 
4> — 4> c , the symmetry breaking field slowly rolls down its 
potential while still inflating the universe, until it starts 
oscillating around oo, ending inflation. The amplitude of 
oscillations is very small compared wi th op, i.e. £ -C 1. 

As we already discussed in Section IV B , in this case 
preheating is inefficient unless there exists a field x 
strongly coupled to cr, which, however, does not acquire 
mass greater than M after spontaneous symmetry break- 
ing. It is especially difficult to arrange it in the case un- 
der consideration, because the amplitude of spontaneous 
symmetry breaking is ~ Mp, which is 16 orders of mag- 
nitude greater than M . Thus one may expect that even 
a very weak interaction of the f ield c r with x wm make 
the field x too heavy. In Section IV B we mentioned that 
it is possible to have light particles which are strongly 
coupled to the field cr, but it requires either fine tuning 
or introduction of flat directions. We will not discuss this 
possibility here. 

Now let us consider a possibility of a perturbative de- 
cay of the oscillations of the field cr. In general, this pro- 
cess could be efficient if this field has a renormalizable 
interaction with some light fields, with a large coupling 
constant a 2 . However, there are two problems associ- 
ated with this possibility. First of all, as we just men- 
tioned, such an interaction typically makes the fields in- 
teracting with the field a extremely heavy, with a mass 
~ oloq ~ aMp. Such particles cannot be produced by 
perturbative decay of the particles with mass M ~ 1 TeV 
unless a is extremely small, a < 10 -16 . But in this case 
the decay rate is extremely small too. 

However, even if one finds a way to have some light 
particles strongly coupled to the field a, its decay still 
occurs very slowly. Indeed, suppose that the decay rate 
r is very large, T ~ M. Then naively one could expect 
that the field decays within the time At ~ ~ M -1 , 
which is of the same order of magnitude as the Hubble 
constant at the end of inflation in this model. In such 
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a case the energy density of the oscillating scalar field 
p a — A/ 2 Cq — M 2 Mp would be rapidly converted to 
the thermal energy ~ T 4 , which would heat the universe 
up to the temperature T r ~ 10~VA/Aip ~ 10 7 'M @. 
However, this is hardly possible. Indeed, strongly inter- 
acting particles in thermal equilibrium typically acquire 
effective mass ~ T. As soon as this mass becomes com- 
parable to M, reheating shuts down. 

As a result, instead of rapidly reaching the temperature 
T r ~ 10 _1 \/-kf Afp, the universe remains at the nearly 
constant temperature T ~ M ~ 1 TeV. The decay of the 
field a continues for a very long time, until the energy 
density of this field drops down from M 2 Mp to AT 4 ; see 
Ref. for a discussion of a similar effect. Thus, reheat- 
ing temperature in this model may be very small; reheat- 
ing may end just before the electroweak phase transition 
or even after it. This may have important implications 
for the theory of baryogenesis and for the cosmological 
moduli problem. 

But the situation may be even more interesting and 
complicated. Indeed, as was noted in Ref. [Q, the ex- 
istence of a second stage of inflation in this model may 
produce a very narrow peak in the spectrum of density 
perturbations. A further investigation of this issue have 
shown that this peak is in fact much higher than origi- 
nally expected, see Ref. [Q. The density perturbations 
corresponding to this peak may lead to a copious black 
hole production even if the universe after inflation rapidly 
enters the radiation dominated stage. However, this ef- 
fect becomes even much more pronounced if the universe 
for a long time remains matter dominated, because with- 
out radiation pressure there is nothing to prevent gravi- 
tational collapse. Therefore even if the peak in spectrum 
of density perturbations is not very high, it may still lead 
to a copious black hole formation. 

If our analysis of reheating in this model is correct, 
then reheating completes at a time which is at least 16 
orders of magnitude longer than the age of the universe 
at the end of inflation. All this time the universe remains 
matter dominated. In this case we expect a huge number 
of black holes with a narrow range of masses to be pro- 
duced soon after inflation, so that most of the matter af- 
ter inflation will be in the form of primordial black holes. 
If the second stage of inflation is prolonged, the masses of 
the black holes will be extremely large, which will lead to 
unacceptable cosmological consequences. However, for a 
wide range of the parameters the second stage of infla- 
tion will be relatively short. In such a case the black hole 
masses will be very small, and they will evaporate before 
primordial nucleosynthesis. The peak in the spectrum of 
density perturbations in this model is very narrow. As a 
result, all black holes will have mass of the same order of 
magnitude. Evaporation of these black holes will happen 
almost simultaneously, at the time related to their mass. 
This process will produce all the radiation and matter we 
see in the universe today 14 1 . 



VI. CONCLUSIONS 

In this paper we considered the initial stages of reheat- 
ing in the hybrid inflation scenario. We have found that 
in certain cases this process begins with a stage of para- 
metric resonance. In such a situation the standard per- 
turbative approach to the theory of reheating should be 
considerably modified. But even in the situations where 
preheating is inefficient, a detailed investigation of the 
behavior of the coupled system of two fields 4> and a is 
necessary for a proper investigation of the perturbative 
reheating regime. For example, as we have shown, for 
A ^> g 2 all energy after inflation is stored in the oscilla- 
tions of the field 4>, so one should study perturbative or 
nonperturbative decay of this field. For A <C g 2 all energy 
is stored in the oscillations of the field a, so one should 
study its decay. Meanwhile for A ~ g 2 the fields enter 
a regime of chaotic oscillations where energy density is 
transferred back and force between the fields <p an d <J - 
Investigation of this regime may be interesting not only 
for the investigation of reheating, but also for the study 
of symmetry breaking pattern in realistic models of hy- 
brid inflation. We have found that preheating can be 
very efficient if the effective masses of the fields cf> and a 
are much greater than the Hubble constant at the end 
of inflation, or if these fields are coupled to other light 
scalar (or vector) fields \- 

In addition to the simplest hybrid inflation scenario |j) , 
we also studied a recently proposed scenario with two 
stages of inflation Jl3|,[l4| . Rather unexpectedly, we found 
that preheating, as well as the standard mechanism of 
reheating, in this scenario is extremely inefficient. In 
such a situation, reheating may occur in a very unusual 
way, via formation of primordial black holes and their 
subsequent evaporation. Even though this process may 
seem somewhat exotic, it may be the leading mechanism 
of reheating in inflationary models with a peak in the 
spectrum of density perturbations |l4| ], as well as in a 
much more general class of hybrid inflationary models 
with a blue spectrum of density perturbations. 
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